Strain sensing with sub— micron sized Al— AlOx— Al tunnel junctions 
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We demonstrate a local strain sensing method for nanostructures based on metallic Al tunnel 
junctions with AlOx barriers. The junctions were fabricated on top of a thin silicon nitride mem- 
brane, which was actuated with an AFM tip attached to a stiff cantilever. A large relative change 
in the tunneling resistance in response to the applied strain (gauge factor) was observed, up to a 
value 37. This facilitates local static strain variation measurements down to ~ 10~^. 



Tunnel junctions with aluminium oxide (AlOx) barri- 
ers are widely used in nanoelectronics applications, such 
as single electron transistors (SET)[l|, SQUIDsd, ra- 
diation detectors 's'l , superconducting quantum bits [1], 
thermometers and coolers [5, 6]. However, their use in lo- 
cal mechanical strain and displacement sensing has been 
very limited, one rare example being measurements of the 
mechanical vibrations of a nanotube[7]- In this letter, we 
show that conventional sub-micron sized aluminum tun- 
nel junctions with AlOx barriers can be used for very 
sensitive strain and displacement detection in nanostruc- 
tures. If used as a displacement detector, our scheme 
is much simpler than the standard scanning tunneling 
microscopy techniques, as the sensor is fabricated on the 
mechanical structure without a need to position and con- 
trol an external electrode with a vacuum gap barrier. 

In typical applications, microscale displacement detec- 
tion is based either on capacitive or optical sensing [8l|. 
Both of these methods have limitations in the sub- 
micron scale: the capacitance of a device can be dom- 
inated by parasitic capacitances of the measurement 
setup and optical sensing is mainly limited by diffraction. 
In addition, there are several other methods for strain 
and displacement sensing based on the phenomenon of 
piezoresistivity[§, [l^, [HI. Typically semiconductors or 
thin metal films are used as piezoresistive transducers. In 
addition, a few groups have reported applicability of mag- 
netic tunnel junctions for strain sensing [l^, [13]. How- 
ever, all methods mentioned above have some limitations 
regarding to possible applications: semiconducting films 
are not very suitable for high frequency applications due 
to their high resistivity , metal films have a relatively 
low response to the applied strain[l^. and magnetic tun- 
nel junctions require external magnetic fields. Also, up 
to date the magnetic tunnel junction sensors have been 
relatively large (20 x 20 Atm2)[l||. 

Here, we stress that lithographically fabricated tun- 
nel junction devices have several advantages in strain de- 
tection, such as their small size (dimensions are easily 
made to be <100 nm), existing high-frequency read-out 
schemes ease of integration into a mechanical sys- 

tem and a good response to the applied strain. In addi- 
tion, the measurement scheme is very simple, as only a 
measurement of the change in the tunneling resistance is 
needed. 

Our experimental setup consisted of a thin (30 nm or 



FIG. 1: (Color online) Schematic view of the experimental 
setup. Zoom-in: SEM image of a tunnel junction. Note that 
the plane of the junction is parallel to the membrane. 



65 nm) silicon nitride (SiN) membrane (dimensions 250 
/im X 250 /ini x 65 nm or 120 /im x 150 fim x 30 
nm), onto which tunnel junctions of lateral dimensions 
200 nm x 300 nm were fabricated with conventional elec- 
tron beam lithography and multi-angle vacuum evapora- 
tion techniques 17]. The membrane was actuated by an 
atomic force microscope (AFM) tip and the tunneling 
resistance was measured simultaneously with standard 
lock-in techniques. In addition to the tunneling resis- 
tance, the photodetector signal and the z-piezo move- 
ment, i.e. the vertical displacement of the AFM can- 
tilever were recorded. All measurements presented in this 
paper were performed at ambient conditions. 

An AFM tip works as a local actuator for the mem- 
brane, and therefore provides the possibility to measure 
the in-plane spatial dependence of the strain field. In ad- 
dition, the AFM method provides a high displacement ac- 
curacy, particularly in the vertical direction. Strain was 
applied to the membrane by pressing it slowly down to 
a predetermined depth in 100 seconds (usually 3/im), 
keeping the maximum displacement for 60 seconds, and 
then releasing it back to its original shape in another 
100 seconds. Similar actuation scheme has been recently 
used to study the elastic properties of thin graphene 
sheets [H, [l^. Stiff cantilevers [13] were used in order 
to minimize the bending of the cantilever during actu- 
ation, and thus exerting most of the force to the mem- 
brane. From the measured force vs. displacement (FZ) 
curves the actual vertical displacement of the membrane 
was calculated by substracting the bending of the can- 
tilever from the z~piezo movement, as described in Rcfs. 
[isl . ITgj . A schematic of the experimental setup and an 
SEM image of a typical tunnel junction are shown in Fig. 
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FIG. 2: (Color online) (a) Tunneling resistance Rt normal- 
ized with unstrained (zero displacement) resistance -Ro as a 
function of vertical displacement. Press and release phases 
are shown with black solid line and red open circles, respec- 
tively. Top curve is data for a 65 nm thick membrane with 
tunneling resistance Rt ~ 26 kQ and bottom curve for 30 nm 
thick membrane with Rt = 9 kf2. (b) Calculated gauge fac- 
tors for the 65 nm (black squares) and the 30 nm data (open 
circles) in (a), (c) Relative change in tunneling resistance vs. 
calculated strain |ezz| for the 65 nm data in (b), (d) same for 
the 30 nm data. Red solid lines in (c) and (d) are calculations 
based on Eq. ([1} with 00 and d as free parameters. 



Examples of the measured response of the tunnehng 
resistance Rt vs. displacement (a lateral distance 1.5 
fim away from the location of the tunnel junction) for 
two different samples is shown in Fig. [H^a). We have 
measured several more samples with essentially identical 
response in the normahzed units Rt/ Ro, where Rq is the 
unstrained resistance. Short-circuited samples showed 
no observable response, thus the observed signal is due 
to tunneling, not due to the piezoresistive behavior of 
the Al metal film The top curve is for a 65 nm 

thick membrane with Rt = 26 kfi, while the lower curve 
is data for a 30 nm thick membrane and Rt — 9 kil. 
Black solid lines correspond to the pressing and open 
red circles to the release of the membrane. The origin 
of the reproducible hysteresis seen for the thin 30 nm 
membranes is mechanical, as it was also verified by the 
FZ curves (not shown here). Detailed understanding of 
the hysteresis is currently lacking. 

To gain more quantitative understanding of the strain 
responsivity of a tunnel junction, we discuss here a model 
based on the widely used simple trapezoidal barrier pic- 
ture (Simmons model) for a tunnel junction [21j and an 
isotropic elastic solid. The change in the tunneling resis- 
tance is caused by a deformation of the tunneling barrier, 
which manifests itself as a change of the effective bar- 
rier thickness d, the barrier height (jjQ and the area A of 
the junction. From geometrical arguments one can easily 



write how the strained values of d and A depend on the 
compressive (diagonal) strain components e{xx,yy,zz} in 
our geometry of a parallel plate junction in the x-y plane 
(plane of the membrane) with the tunneling current along 
the z-coordinate: d = d(l + tzz), A = A{1 + e^x + ^yy)- 
The shear strain components will not have an effect in 
the lowest order. Thus, the strained tunneling resistance 
Rt at zero bias voltage is 



Rt — 



1 



rAi 



■ exp 



y\/2m0o(l 



(1) 

where 4>o is the average barrier height in the strained 
configuration, changing because of the variation of the 
band gap with strain (deformation potential theory) [1^ . 
At first sight the dependence on the plain strains e^xx.yy} 
seems weak; However, from the elasticity theory for 
isotropic solids we can derive that ezz — Tzz{l — v — 
2i^^)/[(l — i')E] — [cxx + ^yy)^l{^ — v), where v is the 
Poisson ratio, E the Young's modulus and Tzz the stress 
perpendicular to the membrane. This means that even in 
the normal-stress-free regions T^z = of the membrane 
(far from the tip position) Rt still has responsivity due 
to the in-plane strains (effect of the finite Poisson ratio) . 

One important figure-of-merit in strain sensors is the 
gauge factor 7, which is defined as a ratio of the rela- 
tive change in the resistance to the applied strain e, i.e. 
7 = (Ai?/i?o)/e- In our experiment, the direct measure- 
ment of the local strain at the tunnel junction location 
is impossible, as no nanoscale commercial strain gauges 
exist. However, we have calculated the strain fields in 
our experimental geometry (SiN membrane and Al wire) 
using the finite element method (FEM) with the displace- 
ment at the tip location as the boundary condition [23| . 
All simulations were performed using non-linear strain 
(large deflection theory), as the displacements in the ex- 
periments are much larger than the membrane thickness. 

As we can see from Eq. ([1]), the tunnel junction re- 
sponds mainly through tzz, so that in our case the most 
reasonable way to define a gauge factor is with that strain 
component. The calculated gauge factors for a number 
of membrane displacements are shown in Fig. [D^b) , using 
the experimentally determined resistance response from 
the data in Fig. [H^a). The 65 nm membrane shows a 
peak value of 7 = 37 at a displacement of 1.5 /.tm, and 
a slightly declining trend with increasing displacement, 
whereas the 30 nm values are clearly lower with a max- 
imum of 7 = 19, but with a different, decreasing trend 
with the displacement. The values of 7 are about an or- 
der of magnitude hig her than for typical thin metal films 
such as gold [l^, [ij . The overall difference between the 
7 values of the two membranes can be explained by the 
difference in Rt, which results from differences in d and 

With the simulated strain, we can also study how the 
measured tunneling resistance varies with e^z, instead 
of just displacement. Figures [5] (c) (65 nm membrane) 
and (d) (30 nm membrane) show this for the same strain 



3 




-120 -80 -40 40 80 120 

Distance (fim) 

FIG. 3: (Color online) Calculated strain \ezz\ (open cir- 
cles) and measured relative change in the tunneling resistance 
(black squares) as a function of distance from the center of 
the membrane through the diagonal of the membrane (white 
dashed line in the inset). Inset shows the calculated strain 
field, when membrane is actuated in the middle. 



simulations as were used in in Fig. [D^b). We see that 
the response is fairly linear in both cases, and that the 
value of strain is much larger for the thinner membrane. 
The quasi-linearity of the response may seem surprising, 
based on Eq. ([T]). We point out, however, that as the ab- 
solute values of strain are ~ 10~^, the exponent is small 
enough to behave linearly. This is verified by plotting 
theoretical points based on the computed strain and the 
modelled resistance response from Eq. ((T)) in Figs. (c) 
and (d), also. Moreover, the different dependence of the 
gauge factor on the displacement can also be explained 
by our model, Eq. In all the theoretical response 

calculations, we estimated the effect of 00 by using a 
typical value for the deformation potential constant in 
semiconductors a = —4 eV, where the change of energy 
gap A_Bg = a{exx + ^yy + ^zz) l23]- The accurate value 
for AlOx is not known to us. 



In addition to local measurements, we have investi- 
gated the spatial dependence of the response by actuating 
the membrane at different points in the plane {x-y coor- 
dinates) of the membrane, and measuring the response 
at a fixed location. The results for the experimental re- 
sponse and for calculated t^z are shown in Fig. [3Ja), 
where the actuation point was moved along the diagonal 
of the membrane (dashed white line in Fig. Oinset). The 
inset shows the shape of the simulated strain field when 
the center of the membrane is pressed. Correspondence 
between the strain field and the tunnel junction response 
is good at a qualitative level, although some differences 
remain that we do not yet fully understand. 

In summary, we have shown that sub-micron sized Al- 
AlOx-Al tunnel junctions can be used for strain sensing 
in low-dimensional nanostructures, yielding relatively 
high gauge factors. By measuring the resistance more ac- 
curately with a realistic sensitivity AR/R ~ 10~^, static 
strains down to 10~^ could be measured with 7 ~ 10. 
That would correspond to a static displacement sensi- 
tivity of 0.4 pm for a cantilever of thickness 1 /im and 
length 1 /im at the cantilever end. We also showed that 
a simple model of the tunnel barrier can explain the re- 
sponse behavior well, and can be used to predict that a 
higher barrier would increase the responsivity. For typ- 
ical AlOx tunnel barriers, 0o ~ 1 eV 2J], which is well 
below the bulk value closer to 10 eV. Thus if one could 
find a way to make barriers closer to the bulk values, the 
responsivity would increase by a factor three. Low tem- 
perature measurements would also help in reducing the 
noise generated by the junction, and would provide more 
information on the microscopic behavior of AlOx barriers 
under strain, because information of the chargin g en ergy 
(capacitance) variation could then be measured |24l | . 
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